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We study four distinct famílies of Gibbs canonical distributions defined on the Standard complex, 
quaternionic, real and classical/nonquantum two-level systems. The structure function or density of 
states for any system is a simple power (1, 3, or -1) of the length of its polarization vector, while 
the magnitude of the energy of the system, in all four cases, is the negative of the logarithm of the 
determinant of the corresponding 2x2 density matrix. Functional relationships — proportional to 
ratios of gamma functions — are found between the average polarizations with respect to the Gibbs 
£ — , distributions and the effective polarization temperature parameters. In the Standard complex case, 

^\ • this yields an interesting alternative, meeting certain probabilistic requirements recently set forth 

Q\ | by Lavenda, to the more conventional (hyperbolic tangent) Brillouin function of paramagnetism 

(which, Lavenda argues, fails to meet such specifications) . 
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PACS Numbers 05.30. Ch, 03.65. Bz, 05.70.-a, 75.20.-g 



I. INTRODUCTION 



In a number of forcefully-written papers some twenty years ago, Band and Park [[Lfj4| strongly recommended that the 
çC) ', bàsic objective of quantum statistical thermodynamics should be taken to be the estimation of a Gibbs distribution 
t-H ' (satisfying imposed constraints on expectation vàlues of observables) over the continuum or "logical spectrum" of 
possible density matrices describing the system in question. This was viewed as a conceptually preferable alternative 
to that of estimating such a distribution over simply a set of eigenstates — corresponding to a single canonical density 
matrix — as in the Standard Jaynesian approach They, however, encountered "two barriers" in attempting 

■ to establish that their suggested methodology would yield identically the same expectation vàlues as the empirically 
successful Jaynesian strategy. "One is essentially philosophical: The prior distribution needed in this continuum 
problem because of the inadequacy of the Laplacian rule of indifference, remains unknown. The other obstacle is 
mathematical: . . . we could not perform the required integrals because we have no uscful analytical description of the 
domain [of density matrices] and its boundary" j|[ p. 235]. More recently, Park (cf. ||), in studying Schròdinger's 
probability relations, wrote that the "details of quantum thermodynamics are presently unknown" and "perhaps there 
is more to the concept of thermodynamic equilibrium than can be captured in the canonical density operator itself." 



The first and principal model considered here (sec HA) is a particular implementation for the case of two-level 



(standard complex) quantum systems of the general program (formulated in terms of systems of arbitrary dimen- 
sionality) of Band and Park (cf. |^]). It yields seemingly novel analyses of both the thermodynamic properties of 
a plane radiation field JÏÏJ and of paramagnetic phenomena (cf. |ÏÏ|- [T^JÏ9[| ) . A noteworthy aspect of our analyses 
is that we find it necessary to identify the parameter /? of the Gibbs distributions introduced below, not with the 
inverse thermodynamic temperature, as is conventional, but rather with the effective polarization temperature JÏ0| , p2| . 
Asymptotically, we find a negative log-linear relationship (equation ( fh?] ) and Fig. Q) between this temperature and 
the reduced thermodynamic temperature of the Standard model of paramagnetism. 

The usual textbook treatment of the particular (spin-1/2) case of paramagnetism in which an ensemble of n 
noninteracting two-level systems is subjected to a magnètic field H predicts that the equilibrium magnetization (M) 
is given by the corresponding (hyperbolic tangent) Brillouin function, [ p"4| , p. 192], 

Motanh^, (1) 

where Mq is the saturation magnetization, /i is the Bohr magneton, k is Boltzmann's constant and T is temperature. 
Brillouin (generalized Langevin) functions, in general, it appears, yield good predictions except for low temperatures 
and high fields |ïq ]. In volume I of his text, Balian |ï^, sec. 1.4.3] writes that the hyperbolic tangent Brillouin "model 
gives a qualitative understanding of the saturation effect... However, the quantitative agreement is not good as the 
behaviour predicted...altogether does not have the same shape as the experimental curves" (cf. Fig. 6.5]). 

Lavenda p. 193] has recently argued that the hyperbolic tangent "Brillouin function has to coincide with the 
first moment of the distribution [for a two-level system having probabilities eX ^ e _ x and e J\_ e - x ], and this means that 



1 



the generating function is Z(x) = coshx [where x = í£y]. Now, it will be appreciated that this function cannot be 
written as a definite integral, such as 

m - \ f = S J^l = Q^I 1/2 (P), (2) 



I v [x) = J?ff V I e ±xt sm v ~ 1 ' 2 tdt (3) 



because the integral form for the hyperbolic Bessel function, 

(s/2)" 

v^r(i + u) 

exists only for v > \. This means that J_i/2(x) = (2/irx) 1 ^ 2 coshx cannot be expressed in the above integral form. 
Since the generating function cannot be derived as the Laplace transform of a prior probability density, it casts 
serious doubts on the probabilistic foundations of the Brillouin function. In other words, any putative expression for 
the generating function must be compatible with the underlying probabilistic structure; that is, it must be able to be 
represented as the Laplace transform of a prior probability density." 

Additionally, as the concluding paragraph of his recent book, Lavenda writes Jh|, p. 198], "Even in this simple case 
of the Langevin function, 

M = M {coth(^) - (4) 

we have witnessed a transition from a statistics dictated by the central-limit theorem, at weak-fields, to one governed 
by extreme-value distributions, at strong-fields. Such richness is not possessed by the Brillouin function, for although 
it is almost identical to the Langevin function in the weak-field limit, the Brillouin function becomes independent 
of the field in the strong-field limit. In the latter limit, it would imply complete saturation which does not lead to 
any probability distribution. This is yet another inadequacy of modeling ferromagnetism by a Brillouin function, in 
the mean field approximation. And it is another illustration of our recurring theme that physical phenomena are 
always connected to probability distributions. If the latter does [sic] not exist, the former is [sic] almost sure 
to be illusory." 

We derive here Standard complex, quaternionic, real and classical/nonquantum counterparts — (p6[), (|52|), ( |56| ) and 
( |67| ) — to the hyperbolic tangent Brillouin function (Q), which are, evidently, not subject to the objections to it that 
Lavenda has raised. The question of whether or not they are superior in explaining physical phenomena, warrants 
further investigation. 



II. ANALYSES OF FOUR TYPES OF TWO-LEVEL SYSTEMS 
A. The Standard complex case 

1. Gibbs canonical distributions 

We consider the two-level quantum systems, describable (in the Standard complex case) by 2 x 2 density matrices 
of the form, 

p=U 1 + z X - ÍV )> G- 2 + ï/ 2 + - 2 <i)· (5) 

2 \x + ly 1 — z J 

We assign to such a system an energy (E) equal (in appropriate units) to the negative of the logarithm of the 
determinant of /?, that is (transforming to spherical coordinates — r,9,(j> — so that x — r cos çí>sin#, y — rsin<^sin# 
and z — rcos8), 

E = -log(l - r 2 ) (6) 

(cf. |ï^]). The radial coordinate (r) represents the degree of purity (the length of the polarization vector) of the 
associated two-level system. The pure states — r = 1 — are assigned E = oo, while the fully mixed state — r = 
— receives E = 0. We, then, have (||), 
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r = ±yl - e~ E . 



(7) 



We employ the positive branch of this solution as the density of states or structure function, 



fl(E) = \J\ - e~ E (0<E< oc). (8) 
(For small vàlues of E, this is approximately y/E.) The integrated density of states is, then, 

/•Eq 

N(E ) = / n(E) = 2(tanlT 1 Ü(E ) - H(E )). 
Jo 

Inverting this relationship (cf. ([ï])), we have 



(9) 



which for large Eq, yields 



^o)=tanh^ + 2 ^ (10) 



n{Eo )^tanh?^±±í. (11) 



The Gibbs distributions having Çl(E) as their structure function are 

-j3E 



f(E;/3) = ^--Cl(E) (0</?<oc), (12) 
where (3 is the effective polarization temperature parameter fïofl and 

z <« = 3ïf^) < 13 > 

serves as the partition function. In general, the zeros of a partition function are of special interest In this 

regard, let us note that the numerator of Z((3) can not assume the value zero, since the gamma function has no 
zeros in the complex plane, although either its real or imaginary parts can vanish ]ïq] . However, the denominator 
of Z((3) assumes infinite magnitude at the [isolated] points f3 = — 3/2,— 5/2,— 7/2,..., since T(x) has simple poles at 
x = —n, n = 0, 1, 2, . . . with residues (— l)"/n! 

We can compute the average or expected value of E, as 

(E) = -| log Z(fl = *(3/2 + 0) - W) = f (14) 

where the digamma function, i/K-z) — Zjkog-'-X 2 ) = r(z/ • (The last equality in (|Ï4|), employing the Pochhammer 
symbol (A)„ is taken from pï| , where it is also shown that this series can be written in terms of a generalized 
hypergeometric 3-F2 series.) For small vàlues of f3, we have 

"" w-a-W (15) 

For large vàlues of (3, 

^WY (16) 

This can be understood from the fact that the digamma function satisfies the functional equation, 

^(l+/?)-^(/?) = I, (17) 

or, more formally, from the asymptotic (f3 — > 00) expansion of (E) — -0(3/2 + (3) — ip(f3), 
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— — ' — — _L n(— 

2/3 8/3 2 + 4/3 3 64/3 4 + 16/5 5 + ( /3 6 



(18) 



(It is highly interesting to note triat one obtains the maximum-likelihood estimator (3 — 3/2(E) for an ideal monatomic 
gas, for which the logarithm of the partition function is — | log/3 |24|, p. 468].) The variance of E — that is, {(E—(E)) 2 } 
— is given by 



var(£) = log Z(/3) = j/tf) - ^(3/2 + /9)« ^ (0 - oc). 



(19) 



Interchanging the roles of E and f3 |23J-|2^], we obtain a (modal, as opposed to maximum likelihood |23|) estimate 
of the effective polarization temperature parameter /?, 



_d_ 

dE 



1 



= 2(e *-l) 



(20) 



Also, we have that 



dE 2 



logíï(E) 



2(e E - l) 2 



(21) 



,2. Approximate duals of the Gibbs canonical distributions 

The square root of ( |Ï9|) is the indicated (minimally informative Bayesian) prior over the parameter /3 |p4]| . It is most 
interesting to note that as [3 — > oo, this is proportional to 1//?. The prior has arisen in several thermodynamic 
models and "is none other than Jeffreys' improper prior for f3 based on the invariance property that the prior be 
invariant with respect to powers of /3" or "the result obtained by Jeffreys for a scale parameter" püfl . The dual 
distribution (now regarding E as the parameter to be fixed, rather than /?) can then be written approximately (due 
to the use of ©) as || 

\z V) — ■ < 22 > 

As an example of this (approximate) duality, we have set (E) = 16.3 in (|2^) and divided the result by .984296 to obtain 
a probability distribution over fi € [0, oo]. The expected value of (5 was, then, found to be .0636579. Substituting this 
value into (|Ï4|), we obtain an (E) of 16.2805 — which is near to 16.3. (For an initial choice of (E) < 16.3, we would 
expect the duality to be more exact in nature.) Now, in general, 

<» - m log zi m » = i<i? w( l + m > - * ( é> » " é> (23) 



and 



d 2 3 
var(/3) « \ ogZ (——) 

+ ^r) - ^)) + 3(^(| + ^r) - V'(^))) « -L «£> - 0). (24) 



4(£) 4V ' ,v,v 2 2{E)' ^ X 2{E)" V1V 2 2(£) y T "2{E)"> 2(E) 

3N 
2<B) 2 



(The variance of /? for an ideal gas of N partides is §| p. 207].) The square root of var(/3), which is 



approximately proportional to j^j (again, in conformity to Jeffreys' rule for a scale parameter), then, serves as the 



prior over (E) 24, eq. 21]. Following [g4], eq. 27a], we can alternatively attempt to express the prior over (E) as the 
product of Q((E)) and the reciprocal of Z( j^e))- We have plotted this result along with jj^ and both curves display 
quitc similar monotonically decreasing behavior. 
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3. Thermodynamic interpretation of information-theoretic results of Krattenthaler and Slater 



We have been led to advance the model (p|)-([l3D on the basis of results reported in |27j and subsequent related 
analyses. In a one-parameter family (denoted q(u)) of probability distributions over the three-dimensional convex 
set (the "Bloch sphere" p8fl , that is, the unit ball in three-space) of two-level quantum systems ([|) was studied. It 
took the form 

»(«) = ,3/ 2r(1 _ M)(1 _ r2)M (-oo<u< 1). (25) 

We note that the Gibbs distributions ( |Ï2| ) can be obtained from ( p5j ) through the pair of changes-of- variables u = 1—/3 
and r — y/1 — e~ E = Q(E) (which is the positive branch of ([?])), accompanied by integration over the two angular 
coordinates (0<é><7r,0<</>< 27r). (We will in sec. [II B| use this pair of transformations also in the quaternionic, 
real and classical counterparts of the Standard complex analysis.) In Cartesian coordinates, q(u) takes the form 

W 2 "») (26) 



7r 3 / 2 r(l - u)(l - x 2 - y 2 - z 2 ) u 

Making the substitutions X — x 2 , Y = y 2 and Z = z 2 , this becomes a Dirichlet distribution 

r(5/2- u) 

7r 3 / 2 r(l - u)y/XVYVZ(l - X - Y - Z) u ' 



(27) 



over the three-dimensional probability simplex spanned by the possible vàlues of X, Y and Z. ("In Bayesian statistics 
the Dirichlet distribution is known as the conjugate distribution of the multinomial distribution... In this sense the 
discrete coherent states are the quantum conjugate states of the Bose-Einstein symmetric number states" |3^, Remark 
2.3.3].) 



The initial motivation for studying the family (25) in |27|| was that it contained as a specific member — <?(-5) — a 
probability distribution (the "quantum Jeffreys prior") proportional to the volume element of the Bures mètric |3~l·| l 
on the Bloch sphere. (It now appears that all those vàlues of u lying between .5 and 1.5, or equivalently (3 E [-.5, .5], 



correspond to such "monotone" mètrics 32 33J . The desirability of monotone mètrics stems from the finding that an 
"infinitesimal statistical distance has to be monotone under stochastic mappings" p. 486].) We were interested 
in p^j (cf. p5| ) in the possibility (in analogy to certain classical/nonquantum results of Clarke and Barron |36|]) that 
the distribution q(.5) might possess certain information-theoretic minimax and maximin properties vis-à-vis all other 
probability distributions over the Bloch sphere. However, considerations of analytical tractability led us to examine 
only those distributions in the family q(u). 

In p7f , the probability distributions q(u) were used to average the n-fold tensor produets (® p) over the Bloch 
sphere, thereby obtaining a one-parameter family of 2™ x 2™ averaged matrices, C, n {u). The eigenvalues of Ç n (u) were 
found to be 

_ 1 r(5/2-u)T(2 + n-d-u)T(l + d-u) n 
V " 2»r(5/2 + n/2-u)r(2 + n/2-u)r(l- U )' Ü ~ U ' ' ' ' ' ' L 2 J 1 8j 



(n-2d+l) 2 fn + l\ 
n ' d = (n+ï) \ d (29) 



with respective multiplicities, 



(Paral·lel analyses have also been conducted using the real analogue, q re ai(u), given by formula ([33J), of the Stan- 
dard complex probability distribution q(u). However, it has proved more difficult, in this case, to give a formal 
demonstration of the validity of the derived expressions for the corresponding eigenvalues and eigenvectors.) 

The subspace — spanned by all those eigenvectors associated with the same eigenvalue X nt d — corresponds to those 
explicit spin states sec. 7.5.j] (3^] with d spins either "up" or "down" (and the other n — d spins, of course, the 
reverse). The 2™-dimensional Hilbert space can be decomposed into the direct sum of carrier spaces of irreducible 
representations of SU(2) x S n . The multiplicities ( p9| ) are the dimensions of the corresponding irreps. The d-tl·i space 
consists of the union of ( n ™£à+i) C0 PÍ es °f irreducible representations of SU(2), each of dimension (n — 2d + 1) or, 
alternatively, of (n — 2d + 1) copies of irreps of S n , each of dimension ^-^d+i) • 
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In |27[ |, an explícit formula, in terms of the eigenvalues, (|28j) was obtained for the relative entropy of p with 
respect to Ç n {u), 

n 

-nS(p)-Tr(®p-logÇ n (u)), (30) 

where S(p) = — Trplogp is the von Neumann entropy of p. Then, the asymptotics of this quantity ( |30| ) was found 
asn-* oo. We can now reexpress this result ( p7L eq. 2.36]) in terms of the thermodynamic variables ((3,E), rather 
than r and u, as (using the positive branch of (pp) 

l logn - \ - \ log2 + PE + ^ log[^§||] + logIX/3) - lo g r(§ + + O(i). (31) 

Let us naively (that is, ignoring the error term immediately above) attempt to find an asymptotically stationary point 
of the relative entropy ( p0|) by setting the derivatives of (|3^) with respect to (3 and E to zero, and then numerically 
solving the resultant pair of nonlinear simultaneous equations. The derivative with respect to (3 simply yields the 
condition ([Ü]) — that is, the likelihood equation — that E = (E) , while that with respect to E can be expressed as 

'^^^^m'^TiW)^ (32) 



We obtained a stationary value of the truncated asymptotics (J3Ï]) at the point {(3 — .457407, E — 2.58527). It appears 
possible that this point serves as an asymptotic minimax for the relative entropy. Though additional numerical 
evidence appears consistent with such a proposition, a formal demonstration remains to be given. In pï| , the 
asymptotics of the maximin — with respect to the one-parameter family — of the relative entropy was formally 
found. (This involved integrating out the radial coordinate r. In this analysis, it was possible to show that it was 
harmless to ignore the error term in (|3Ï|).) It corresponded to the value (3 = .468733. This was obtained as the 
solution of an equation |27], eq. (3.10)], here expressible as 

2f3 3 v&r(E) = 1. (33) 



(Use of var(_E) rss from (19) in this equation would give a solution of f3 = |.) It is interesting to note that both 



the [rather proximate] vàlues of (3 of .457407 and .468733 lie in the range [-.5, .5] associated with monotone mètrics 
|j32f . In their classical/nonquantum analysis, Clarke and Barron [|36| found that both the asymptotic minimax and 
maximin were given by the very same probability distribution — the Jeffreys prior, which is based on the unique 
monotone mètric in that classical domain of study. 

Let us note — in view of the forms of expression of (31) and ( j32|) — that in [jï(], eq. (2.4)] (cf. 22|) the expression 



1 A + íl(E), 1 r l + 7\ r 3 r 5 

was taken to be the reciprocal (t _1 ) of "an effective polarization temperature (r) [which] should not be confused 
with the radiance temperature obtained using Planck's spectral law" . (The equality in |I1) will not hold in the 



quatcrnionic and real cases discussed below (sec. II B), since then the structure function will not simply equal the 
radial coordinate.) Inverting equation (2.4) of |Ï0[] , we obtain 



e 2 / T - 1 1 

= tanh — (35) 



e 



(cf. formulas (1-28) and (1.36) of Q). It will be of substantial interest (setting r = f3) to compare ( p5|) with the 
relations we obtain — (|3^) , ( |52| ) and (||^) — for the average polarization (r) as a function of (3. This will be done in 
Fig. |. 



4-. Dependence of the average polarization upon the parameter j3 

The expected value of the structure function Q(E), that is (||) — equal to the radial coordinate (r) — with respect 
to the family of Gibbs distributions (113) is exactly computable. We have that 
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(36) 



For small vàlues of (cf. (U 



Asymptotically, —> aa, 



1 2 5 73 575 1 

(r) ~ 0F ( /?V2 ~ 4/?3/ 2 + 64/^T 2 " ~ 5120V 2 ' + 



We have found a series of curves for integral n that converges (as C. Krattcnthalcr has formally demonstrated) to 
( |36| ) as n ~ * oo (cf. |39]]). The curves are obtained by summing the weighted terms (n — 2d)/n over the 1 + |_n/2j 
subspaces (indexed by d) of explicit spin states, using as the weights, the probabilities assigned to the subspaces, that 
is, the product of the eigenvalue (\ n ,d) for the d-th subspace (|8|) and its associated multiplicity (m n ,d) (p9[). Thus, 
we consider 

L"/2J _ , 

( )TOn,dA n ,d· (39) 

n 

This gives the average or expected polarization (since n — 2d is the net spin for n — d majority spins assigned +1 and 
d minority spins assigned —1), ranging between and 1. (There have been recent studies f4Ï| , [42|| of the temperature 
dependence of the average spin polarization in a system — a quantum Hall ferromagnet — in which the spins are not 
[as assumed here] noninteracting (cf. |^3| )-) Thus, similarly to the parameter t in |ïo[| , the value = corresponds to 
complete polarization (total domination by the majority spin) and = oo to the case in which there is no predominanec 
of spin in any direction. 



The convergence of the sum (39) to the mean value (B6J) has been essentially demonstrated by Krattenthaler by, 



first, rewriting ( p9| ) as the difference 

h/2J n _ 2d 1 L»/2J 1 

V" ( )m n , d \ nid - (-)m n d X n d . (40) 

* — ' n í — ' n 

d=Q d=0 

(His proof — applicable to any probability distribution over the Bloch sphere which is spherically-symmctric is 
conducted in terms of the spherical coordinates used in pTj], but it easily carries over to the presentation here in 
terms of the thermodynamic variables, and E.) He shows — employing Theorem 15 of [^7j to reexpress X n ,d, thcn 
splitting the rcsulting expression into four terms, in the manner of (2.41)-(2.46) of p7|] and applying the binomial 
theorem — that ( ^p| ) and, hence, (^) can be rewritten as the sum of the expected value of q(u) (that is, J^r^/^-u) ) 
and a term that converges to zero as n — > oo. 



5. Relations between average polarization ( jg^ ) and hyperbolic tangent Brillouin (0) functions 

Let us, in the obvious task of trying to relate the results here to the conventionally employed hyperbolic tangent 
Brillouin function (|ï|) , equate the reduced magnetization |ïï| given by -j^ — tanh ^ to the average polarization (r) 

obtained from formula ([56]). To do so, we set the reduced temperature [x = ^) equal to tanh _1 (?'). Asymptotically 
then, as (3 — * oo, this gives 



0?%/^ 12tt 3 / 2 /3 3 / 2 y (3 2 
Consequently, for large vàlues of 0, we have the log-linear approximation, 

logx w log2 - l - logTr - l - log/3 = .120782 - i log/3. (42) 

In Fig. ^, we plot logtanh _1 (r) vs. the logarithm of the effective polarization temperature 0. For large 0, the curve 
is approximately linear, having a slope of — \. 



7 



log inv. tanh 




FIG 
of sec. 

linear with a slope of — i. 



The Logarithm of the Inverse Hyperbolic Tangent of the Average Polarization ( |36[ ) for the Standard Complex Model 
[IA| vs. the Logarithm of the Effective Polarization Temperature Parameter /3. For large 0, the curve is approximately 



To further study the relationship of the results in this study to the Standard treatment of paramagnetism (see 
also the formulas we obtain for the integrated density of states — (|ï(]), pq), d65|)), in which the hyperbolic tangent 
Brillouin function ([]]) is employed jÏT might prové useful to employ the relationships between the hyperbolic 

functions and the gamma functions in the complex plane eqs. 6.1.23-6.1.32]. Pursuing this line, we have found 
that 



20F(/3 + l)r(3/2 - (3) 
(402 - l)T(-0) 



= i tanh —r- = tan ptt. 



(43) 



Error bounds for the asymptotic expansion of the ratio of two gamma functions with complex argument (cf. ( p^ ) and 
(pq)) have been given in Eq] (extending earlier work Bfl] on the ratio of two gamma functions with real argument), 
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using generalized Bernoulli polynomials (cf. Mj). 



B. Cases other than the Standard complex 

Although quantum mechanics is usually treated as a theory over the algebraic field of complex numbers, it is also 
possible and of considerable interest to study versions over the real and quaternionic fields, as well Q,^9). The 
possibility of testing these different forms against one another has been studied ]50|-|5^|. It appears that the results 
here — giving distinct Gibbs distributions in the three quantum mechanical scenarios — may present another avenue 
to addressing such fundamental issues. 



1. Quaternionic case 

For the quaternions, the analogue of the probability distribution q(u) given in ( p5|) takes the form (cf. eq. 26]), 

. . r(7/2-u)r 4 sin 3 6>ism 2 6> 2 sin6l3 . fAÀ . 
W(") = ^ii-u^-^Y (-oo<«<1), (44) 

being defined over the /ïue-dimensional unit ball or "quaternionic Bloch sphere." Proceeding with the same transfor- 
mations as in the Standard complex case — that is, u = 1 — and r = vT— e~ E — and integrating over the four 
angular coordinates, we arrive at a family of Gibbs distributions (|Ï2] ) having as their structure function or density of 
states (cf. (§)), 

Ü quat (E) = n 3 complex (E) = (1 - er E )i (45) 

(which is approximately E z / 2 for small E), where fl C om P iex(E) = íl(E) = v 7 ! — e~ E , from (||). The integrated density 
of states (cf. @) is, then, 

N qaat (E Q ) = ü quat (E) = 2(tanh- 1 íl complex (E ) + ~ 4e g °) ^ (4g) 

Also, we have the difference, 

Ncomplex(Eo) — N quat (E ) — complexió) ( _ ^ 

As Eq — ► oo, it approaches from below, the limiting value of |. For large Eq (cf. (|j), Jïï])), 

<~» /-ei \ i. 3N qua t(Eo) + 8 

^compíex(^o) « tanh ^— i . (48) 



The partition function (cf. (13)) in the quaternionic case is 

Z q uat((3) = 4r(5/2 + /3) - (49) 

We now have (cf. (jïj)) 

(Vi) = log = ^(V2 + « - í(^) « ^ (/3->oo) (50) 



var ïíiot (B) = ^log^ uat ( y 9)=V'(^)-V' , (5/2 + í 9)« i (0 °°)- (51) 



and (cf. (|19| 



It, thus, appears that the quaternionic results relating to approximate duality can, in essence, be obtained from the 

5 

2 



(standard complex) ones given above (sec. HA2), simply through the substitution of | for |, as appropriate. We 
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have conducted a test of approximate duality similar to that reported (for the Standard complex case) immediately 
after (p2|). We have set (E) quat = 16.3 in the evident quaternionic analogue of fl22]), and divided the result by .902062 
to obtain a probability distribution over (3 £ [0,oo]. The expected value of j3 was, then, found to be .0664174. 
Substituting this value into (|50|), we obtain an estimate of (E) quat equal to 16.2645 ~ 16.3. 
For the average quaternionic polarization, we have 

/ \-/o - 8r(5/2 + /3) _ 8r(5/2 + /3) 

{ quat) ~ ^ at ^> > - 3/3(1 + /3)(2 + /3)(3 + (3)V^r(p) ~ 3^(3 + /?) ' 1 j 

(It would appear problematical, however, to develop a quaternionic counterpart to the analysis of Krattenthaler and 
Slater (^7), as there are difficulties in defining an appropriate extension of the tensor product |JÏ8| .) 

2. Real case 

We can also, proceeding similarly, take 

, , 2r(2-w)r (l-u)r 

^l{u) = T{l _ u){l > r2)u = W ^ W (-00 <«<!), (53) 

as a probability distribution over the unit (two-dimensional) disk. Integrating over the single angular coordinate and 
using the transformations, u — 1 — (3 and r — \/\ — e~ E , we transform ( [53] ) to a Gibbs distribution, having a structure 
function Q rea i(E) = 1 and a partition function, 

ZreM = = p. (54) 



(The integrated density of states is, of course, then, simply N rea i(Eo) = Eq.) Now (cf. (17)) 



(Ereal) = - ^ log Z real (/?) = ^(1 + 0) - l/>(0) = ~ (55) 

So, it appears that the "fraction" ' 1 = | plays the analogous role in real quantum mechanics as | in the Standard 
complex version and I in the quaternionic form. (These results appear consistent with the equipartition of energy 
theorem |33| in classical statistical mechanics, as two, three and five are the corresponding degrees of freedom, that 
is, the number of variables needed to fully specify a real, Standard complex or quaternionic density matrix.) 

The expected value of r with respect to (|||) (again using u = 1 — /3) gives us the average polarization in the real 
case, 

Vir(i + fl 



(57) 



The three forms of average polarization — ( |36[) , (^2|) and ( |56|) — are particular instances of the formula, 

T(l + m/2)r(l/2 + 13 + m/2) 

r(i + /? + m/2)r(i±2i) ' 

where the quaternionic case corresponds to m — 4, the Standard complex to m — 2 and the real to m = 1. 

3. Classïcal/nonquantum case 

The case m = is, essentially classical/nonquantum in nature, corresponding to the family of probability distribu- 
tions, 

, 2T(3/2-u) . /co . 

«-».(«) = ^r {1 _ u){1 _^ (-»<«< D- (58) 



(If one makes the transformation r = V ^i?, this becomes (cf. (|27|)) a family of beta distributions. Then, the member 
corresponding to u = 5, is the (bimodal) arcsine (alternatively, cosine distribution, corresponding to a process 
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with a single degree of freedom, to which Boltzmann's principle is inapplicable p5[ . The arcsine distribution serves 
as the noninformative Jeffreys' prior for the Bayesian inference of the parameter of a binomial distribution ]56| , |57| ] . ) 
Making the same substitutions as in the previous cases, that is u = 1 — j3 and r = y/1 — e~ E , into (psf), we obtain 

^class (E) = - 1 —— = 1 , (59) 

W0 = j^My, (60) 

and 

(E class ) =-^logZ cíass (/3)=^(l/2 + /3)-^(/3)« Yp O 9 ^ 00 ) ( 61 ) 
(corresponding to the single degree of freedom in this case). Also, setting m = in (|57|), we have 

_ r(i/2 + /3) 

IW! '"V^(1 + J3)' 
For the integrated density of states, we find N c i ass (Eo) = 2tanh _1 Q C ompiex(Ea), so 

^complex (Eo )= tanh — £^£ií — (63) 



C. Comparison of results for the four different cases 

The average polarization indices for the four models considered above all equal 1 for /3 = and monotonically 
decrease to for (3 = oo. However, for [3 strictly between and oo, the quaternionic value is always greater than the 
Standard complex one, which, in turn, is always greater than the real index (which itself dominates the classical one). 
(In this regard, it may be helpful to note that the ratio of ( [Ü2"|) to ( |36|) simplifies to 2 ^2+p) ■ ) 

In Fig. H, we plot the average polarization as a function of (3 for the quaternionic, Standard complex, real and 



classical cases (m = 4, 2, 1 and 0, respectively). We also display the Brosseau-Bicout relation (35), having set r = (i. 
This last curve crosses the other four. It predicts greater average polarization in the vicinity of r = = and less 
above a certain threshold (which is smallest in the quaternionic case and largest in the classical one). 
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beta 



FIG. 2. Average Polarization as a Function of the Effective Polarization Tem per ature (3. The quaternionic curve ( |52| ) domi- 
nates the Standard complex curve which, in turn, dominates the real curve (pq). This dominates the classical/nonquantum 
curve (p2[)). The function of Brosseau and Bicout (pg), setting r = (3, making it equivalent to tanh -|, crosses the four other 
curves, intersecting at /3 = .76007, 1.04585, 1.46249 and 3.1857. 



In Fig. ||, we plot the expected value (E) of the energy E for the quaternionic, Standard complex, real and classical 
cases, while in Fig. ^, we show the corresponding variances var(_E).. 
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FIG. 3. Expected vàlues (E) of the energy E. The order of dominance is {E qua t) > (E cornp i ex ) > (E rea i) > (E c i ass ). 
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FIG. 4. Variances of the energy E. The order of dominance is v&r(E quat ) > var(E comp iex) > vax(E rea i) > vax(E c i ass ). 
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D. An alternative family of Gibbs canonical distributions for the Standard complex case 



In thc analyses reported above, we have relied upon certain transformations of the probability distribution q(u) 
(p5|), defined over the Bloch sphere of two-level systems, and its real ( |53|) and quaternionic ( fii]) anàlogues, to obtain 
the Gibbs distributions we have studied. The distribution q(u) has the attractive feature that for u <E [.5, 1) (and 
also for u € [1, 1.5], although q(u) is improper in this range), it is proportional to the volume element of a monotone 
mètric on the two-level systems. (The value u = .5 corresponds to the minimal [Bures] monotone mètric and u = 1.5 
to the maximal one.) 

Nevertheless, there appear to be alternatives to q(u) that are of certain interest and might conceivably be physically 
meaningful. One such family of probability distributions can be expressed as p7| , eq. (3.11)] 

, . (l-u)r(3/2-u)rlog[(l+r)/(l-r)]sin0 . feA . 
Vkmb(u)='- l 2 W(l-.)(l-r^ (-00 <»<1). (64) 

For u = |, this is proportional to thc volume element of thc Kubo-Mori/Bogoliubov (monotone) mètric |34l , |58| , p9| . 
"This mètric is infinitesimally induced by the (nonsymmetric) relative entropy functional or the von Neumann entropy 
of density matrices. Hence its geometry expresses maximal uncertainty" |M - The family (Q) can essentially be 
obtained from the family q(u) by replacing a factor of r by log[ faí^A 1 an< ^ renormalizing the result. For vàlues of r 
in the neighborhood of 0, q(u) and qKMB(u), for any fixed u, are approximately proportional (cf. (j34|)). (The KMB- 
metric is the extreme element [a = ±1] of a one-parameter (a) family — distinct fro m (p4[ ) — all the members of which 
correspond to monotone Riemannian mètrics induced by quantum alpha-entropies |59(| . However, the corresponding 
volume elements for — 1 < a < 1 do not appear to be exactly normalizable over the Bloch sphere of two-level systems.) 

If we perform the transformations, u = 1 — /3 and r — y/T— e~ E , then (|64| ) becomes a family of Gibbs distributions 
( |Ï2| ) with the partition function, 

Zkmb{P) - - mi/2 + í3 y ( 65 ) 



and, quite interestingly, the structure function (cf. (|34| 



ÏÏkmb(E) = log[ 



(i + VT^ 
(i-vT^ 



e- E ) 



= 2(VÏ~^ + [ \ ' +-..) = 2(ü complex (E) + " U f ] +.-.)• (66) 

Using the positive branch of (Q), ÍIkmb(E) is simply twice the reciprocal (t _1 ) of the effective polarization temper- 
atura (r), as defined by Brosseau and Bicout |ï(], eq. 2.4] (cf. (|34"|)). The range of possible vàlues of this structure 
function is the nonnegative real axis [0, oo], in contrast to the structure functions considered above, which only assume 
vàlues on the unit interval [0,1]. We were able to explicitly compute the corresponding integrated density of states — 
but the expression MATHEMATICA 3.0 yielded was highly complicated. In Fig. ^|, we have plotted the five integrated 
density of states functions obtained here — N comp i ex (E ), N quat (E ), N rea i(E ) = E , N c i ass (E ) and N KMB (E ). 
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FIG. 5. Integrated Densities of States for Five Different Models. The most steeply-rising curve corresponds to Nkmb(Eo), 
while the order of dominance for the other four curves is classical > real > Standard complex > quaternionic. The KMB-cmve 
crosses the classical one at Eo = 1.57565, the real linear one — N rea i(Eo) = -Eo — at .53341, the Standard complex 
curve at .000111286, and the quaternionic curve at .0000405489. Asymptotically {Eo — > oo), the difference (^), that is 
N comp i ex (Eo) — N qU at(Eo), monotonically increases to its limit |. 



1G 



The first term of the asymptotic expansion ((3 — > oo) of — ^ log Zkmb is ^g. This is the same as was found above 
flq) for the Standard complex case and which coincides with the behavior of the ideal monatomic gas. The first term 



of the asymptotic expansion (j3 — ► oo) of the square root of JLp log Zkmb is nfy, again as in the Standard complex 



case considered previously fll9|). 

The expected value of r with respect to ( |64| ) — giving us the average polarization index — can be expressed as 

, . 2/^(1/2 + /3) P F Q [{1, 1,2}, {§,2 + /?},!] 

{rKMB) = ^Fr(2 + /3) ' (67) 

Due to its hypergeometric character, it is difficult to analytically study (|67| ). A plot (Fig. ||) of the difference between 
it and the alternative (standard complex) average polarization index (p6|), shows that the latter index is dominated 
by j67|), with the greatest difference (w .0526) occurring in the vicinity of j3 = .49825. (We have earlier indicated 
that for the distributions q(u), given in (p5|), taking u = 1 — fi, the range j3 G [—-5, .5] — which includes .49825 — 
parameterizes a continuum of monotone mètrics.) 
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FIG. 6. The average polarization given by (tkmb) mimis the average polarization (r) — given by (^6|), plotted in 
for the first/primary Standard complex model considered here (sec. II A). 
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III. CONCLUDING REMARKS 



The Brillouin functions of paramagnetism serve as the basis of the mean-field theory of ferromagnetism |^,|L2j,|20| . 
It would, therefore, be of interest to construct mean-field theories based on the alternative functions that we have 



developed above. In this regard, following the usual line of argument |2fj, sec. 6.2] |60, sec. 9.2], we have replaced (3 
by (3/ (A(r)) in the small-/3 approximation (^7|) to the average polarization index (|36|) in the Standard complex case. (A 
corresponds to the molecular field parameter.) Then, we find that there exists a nonzero critical effective polarization 
temperature, 

A = 4(21og2-l) = - 647175A (68) 

below which spontaneous magnetization is possible. We are (using the same approximation and following J6Ï[ sec. 
3.1]) able to construct an order parameter of the form 

2<r)-l = ±(l-^) 1 /2. (69) 

So, the exponent for the power law behavior of the order parameter is the same as for the hyperbolic tangent Brillouin 
function (|ï|) |jsï| eq. (3.10)], that is the simple fraction | — as is not actually the case for real ferromagnets. (Lavenda 
and Florio ]62[ assert that the probability density centered about the metastable state of the order parameter for the 
mean-field theory of the kinetic Weiss-Ising model is an asymptotic distribution for the smallest value, rather than a 
normal distribution, as generally assumed.) 

The possibility of obtaining results analogous to those derived here, for n-level systems (n > 2) should be investi- 
gated. Such results could, then, be compared with those given by the particular n-level form of the Brillouin function. 
In this regard, it might be of some relevance to note that Brosseau has recently shown that the polarization 
entropy of a stochastic radiation field depends on (n — 1) measures of the degree of polarization of the field. 

If one considers the pure states in a 2m-dimensional complex Hilbert space, endowed with the unitarily invariant 
integration measure pjj , then tracing over m degrees of freedom, one induces a probability distribution on the two- 
level Standard complex quantum systems taking the form of the Gibbs distribution (|I|) with (3 = m — 1 |5ï| (cf. 

!§)■ 

Boltzmann's principle — S = k log + const — relates the entropy, S, to the logarithm of the thermodynamic 
probability, O, of a given state sec. II.4.1]. The question of whether this important principle has a meaningful 
application to the structure functions considered here (Acomplexi Qquat, ^reai = 1, ^dass and ÍIkmb) certainly mèrits 
investigation. (We note that all but the last of these functions assume the value 1 for E = oo, the logarithm of which 
is equaling, as would seem appropriate, the von Neumann entropy, — Trplogp, of a pure state — that is, one for 
which, detp = or, equivalently, r = 1 or E — oo (cf. (^)). However, the von Neumann entropy of the fully mixed 
state, corresponding to r — or E = 0, is log 2, not — oo, as a naive application of the principle would, then, give.) 

White |67| 1 has developed a density matrix formulation for quantum renormalization groups and applied it to 
Heisenberg chains. (He shows that keeping the most probable eigenstates of the block density matrix gives the most 
aceurate representation of the system as a whole, that is, the block plus the rest of the lattice.) We have indicated here 
a different use of the density matrix concept in studying spin systems, in line with the program expounded by Band 
and Park. In their series of papers jï] |j, they take exception to the Jaynesian strategy of maximization of the von 
Neumann entropy subject to constraints on expected vàlues of observables j^J(| and suggest that rather one should 
estimate a Gibbs distribution over the continuum or "logical spectrum" of possible density matrices. We should also 
point out, however, that Band and Park — in contrast to Lavenda |24|,|25| and Tikochinsky and Levine [M — do 
not discuss the notion of an (improper or nonnormalizable) structure function or density of states over the energy 
levels (such as we encounter here with (pj), jjq), (p9|), d6Ü| ) and Q. rea i(E) = 1), but solely that of a reparameterization- 
invariant prior probability distribution over the convex set of quantum systems (cf. ]5ffl ) . Within the framework of their 
analysis, concerned with systems of arbitrary dimensionality (not simply two-dimensional systems, as here), Band 
and Park did conclude ||^| that for the case of "strong" or dynamical equilibrium, an axiom of "data indifference" 
was preferable to one of "state indifference" . 

In conclusion, let us state that we have found counterparts (that is, the expected vàlues (r) of the radial coordinate 
r in the Bloch sphere-type representations of the two-level systems) — d36|), (p2|), ( |56| ) and ( |67| ) — to the hyperbolic 
tangent Brillouin function ([!]), which are evidently not subject to the methodological objections of Lavenda H]. Of 
particular novelty has been the identification of the parameter (3 of the Gibbs canonical distributions introduced above 
with the effective polarization temperature, and not, as is conventional, the inverse thermodynamic temperature. 
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Further theoretical and empirical analyses pertaining to the potential applicability to physical phenomena of the 
results reported here would, indeed, appear appropriate. 
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